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۱ I. INTRODUCTION 

۱۳ or ۷ ۰5 و‎ notior cen play a significant role in 
the mecharical energy budget of the ocean boundary layer. 
Linden (1975) carried out experiments To cbserve mixed layer 
deepenirg in a stratified fluid. He felndo that latessal 
waves are important in the energetics of the mixed layer and 
that the mixing process may provide a significant source of 
energy for the generation of internal waves. In this thesis 
the interior motion is prescribed and nor affected by the 
mixed layer processes. However, the modelling of the inter- 
action of fres internal waves with mixed layer processes is 
ar impcrtant area for future study. Linden points cut that 
he did not acccunt for waves generated elsewhere which 
propagaze into the region of consideration. Alsc, the 
internal waves corsidered by Linden were of relatively high 
frequency. 

De Szoeke (1980) found wind stress curl-driven vertical 
advecticn to be as important as the entrainment rate in 
mixed layer dynamics. Cushman-Roisin (1981) likewise found 
wind-driven advecticn to be as important as wind-driven 
Mixing in mixed layer dvnamics and temperature front 
formation. 


Much work has been done to study the ocean response to 


storns. Price et al. (1978), in their work or the mixed 
layer resronse  *o stcrms, found it necessary to include 
vertical advection. Price (1981) found that upwelling wiil 


Significantly enhance the sea-surfase temperature response 


n 


to a slowly moving hurricane. Adamac 3* al. (1981), using 
two-dimensional model to study the cceanic response to a 
hurzicane, found a strong interdependence between mixing and 


advection. sreatbatch (1983) also found advecticn +o be 


11 





important in increasing maximum cooling due to the passage 
@f a storr. 

Only recently has much been done to include long period 
interior motion in mixed layer models. Stzvenso 
(1980,1981) =xanineä the response of the mixed la 
02 ۲+ ۱1۳660 by linear Rossby waves. He found no mixed lays 
response during the first half of the heating season while 
the mixed layer is shallowing and a maximum sea surface 
temperature response at the end of the heating season and 
beginning of the cooling season. Stévenson found that the 
maximum response of mixed layer depth to the wave occurred 
at or near the end of the cooling season. Also, he found 
that the lowest frequency waves produced the greatest mixed 
layer depth and temperature responses. This study differs 
from that of Stevenscn in that a finite amplitude prescribed 
interior motion is used rather than a lirear Rossby wave. 
This produces a phase dependen*+ asymmetry in che mixec layer 


response to the interior motion. 


Burger (1982) . successfully introduced tidal period 
vertical advection into the Garwood (1976,1977) one- 
dimensional mixed layer model. He found that the terior 


in 
motion had a significant effect on the near-surface tempera- 
ture field and that the inclusion of the interior motion can 
improve tke model for single station forecasting. 

In this study, the Garwocd (1976,1977) one-dimensional 
ccean planetary boundary layer model is used. Vertical 
advecticn was added tc the model. The annual heat cycle is 
a sine function and wind stress is held constant at one dyne 
per square centimeter. By keeping the surface boundary 
conditions simples, the effect cf the interior motior can be 
nore readily observed. 

The interior Vertical motion is a prescribed, forced 
interior wave which is linear with depth and sinusoidel in 


time: 
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where D is a reference depth, 200 meters, and Wo is the 
amplitude ct the prescribed vertica motion at depth, D. 
The periods, T, studied are one y2ar, one-half year, onre- 
third year and one-quarter year. The phase difference, $., 
between the interior motion and the annual heating cycle is 
۱ ۱:0 from 0° to 360º at increments of 22.5”. This 
interior motion is general and may ba representative of a 
number cf possible long period phéencmena, for example, 
internal Rcssby waves or perhaps Ekman pumping on seasonal 
scales. 

The oceanic variables of interest are mixed layer depth, 
h, sea surface temperature which equals mixed layer tenpera- 
ture, Tn, the heat content, H, ani the potential energy, P. 
These variables are contoured as functions of time and phase 
fcr each cf the pericds studied. Also, contours aze made of 
the differences of these variables relative to the case of 
no intericr moticn. These contours are used to observe the 
effect of the interior motion on the mixed layer processes. 

Fur-hermcre, the phas2, 9, in equation (1.1) can be 


replaced ty the expression: 


= 9 (1.2) 





where x is horizontal displacement and L is some typical 
wavelength. Then, at any fixed x, squation (1.2) represents 
the vertical motion due to the passage cf a travelling wave. 
The wavelength, L, might be, for example, 300 km for an 
annual karoclinic Rossby wave in thé central north Pacific 
(Kang and Magaard, 1980). The phase in ¢quation (1.1) can 
thus ke related to herizontal location. TEE CONCOURS Of h, 


Tair H anā P versus phase and tims then indicate horizontal 
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A. THE MODEL 


1. Some General Considerations 
Ecr *his study the Gazwood  (1976,1977) one- 
dimensional bulk mixed layer model is used. The model is 


q 


based on the Navier-Stokes equations of mction, assumes the 
continuity ¢quation for an incompressible fluid and makes 
the ons-dimensicnal approximation. The Garwccê model 
differs from prsvious bulk mixed layer models in several 
features. The time scale for viscous dissiparicn cf turbu- 
len= kinetic energy, TKE, is dependent on planetary rota- 
eon. This allews for 3 cyclical steady state on an annual 
basis by increasing dissipation for desper mixed layers. 
Both downward and upward turbulent heat fluxes are important 
RR this model. Furthermore, entrainment depends on the 
relative distribution of IKE between vertical and horizontal 
components. 


Befcre continuing, a comment on notation is in 


order. An overbar cn 2 variable will indicate the horizon- 
tally averaged mean. A prime on 2 variables will indicate 
the fluctuation aboue this @herizental mean. Angled 


brackets, < >, will indicate vertically averaged, mixed 


layer mean. 
2. Eguation of Siate 


In the Garwood model, “he squation of state has the 


fOr: 
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26 
p=p[1-«(T-T,) + 6(s-5,)]. — 


The thermal expansion coefficient, %, éd dore ion OF T 
end is re=computed daily. The constants a, T, and S, are 
typical density, temperature and salinity values,  respec- 
tively. In this study, salinity is fixed (precipitation and 
evaporaticn are disallowed) and sat 2qual to the value of 
S,. Equation (2.1) then becomes 


p=psf1-4(T-T,)| l (2.2) 


B. THE HEAT EQUATION 


The time rate cf change of temperature is given by 


the one-dimensional heat equation: 


de e Q ۳ ee 2 F 
cs _ (re (E) mz «T. (23) 
9 dz د‎ NE dt 22 


Mae cirst term on the right is vertical curbulen+ heat flux, 
the second is solar insolation, che hizo às vertical 
temperature advection and the last, diffusion. 
solar insolation will be assumed to bə abscrbed at the 
surface (in the first few centimeters) of th= ocean. Hence, 
the sclar insolation term will be dropped in favor of a 
ee boundary .cenaitiz on the turbulert heat flux. 
Equation (2.3) then becomes 





3T (2.4) 


a 
| 
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ct‏ 
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In this model, mixed layer temperature is assumed 


to be vertically homcgeneous. Equation (2.4) then reduces 
£0 
7 \ 
CALO " 2 ( L2 (2.5) 
ابت‎ w | 

dt ± 
for mixed layer :emperature, T,,. Obviously, Ty, is equal 
to sea surface temperature because of the vertical homoge- 


neity assumption in the mixed layer. 
The surface bcundary condition on heat flux is given 


by: 





CUm) = 一 a (2.6) 


where the Zer) subscript indicates the value of w’T’' at z=0. 
The net surface heat flux, Q/ P. 2p» is prescribed to be 


hermonic in times 


a - Lo Sin (>= ) , (2:7) 
م6‎ Ac p 








The amplitude of this sine function is set equal to a 
representative mid-latitude value If 4+1073°C cm/séc. In 
this study, heat flux is positive downward. Time t=0 is the 
start of the "heating season", when Q first becomes 
positive. 


The boundary condition on heat flux at 


th ct 

ee! > 
0 
o 
لم‎ 
un 
0 
O 
th 


the mixed layer is given by th? entrainment heat 
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w ! xr ), nd Le ANE bys J 
Ee “ho sobscript, h, indicates the value of w’T’ at 2=-h, 
the basa of “he mixed layer. The variable, We i 


* 


entrainment rate and AT is the temperature jump 


ac 
entrainment zone at the base of thə nixed layer. The turbu- 
lent heat flux is linear with depth in the nixed lay2r (due 
to the assumption of a well-mixed layer), monotonically 
Gecreases with depth in the entrainment zone and is Zero 
below. Ths wn ir 


0 
Foot turbul=r: heat flux is sho 
u 


P 
1 


he net surface hat fl 
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Figure 2.1 Turbulent Heat Flux Profile. 


to negative values according to equation (2.7). The 
entrainment heat flux term, -AT de, is always nega:iv= or 
zer). This is because the entrainment rate is positive when 
the mixed layer is deepening buz "turns oft" and is zerc for 
non-entraining shallcwing cases. ۳ لا‎ ekness otf ë -he 
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entrainmert zone, d, ¿s assumed zo be mucn less *ha- -ns5 
depth of the mixed layer, h. 
Equation (2.5) now becomes 











DNI ل‎ G eve (2.9) 
gt Po Ca N h 
۱ 
Nota that the advecticn term is absent. PAUSE, .hnowanve-cior 
motion does not directly affect tha mixed layer temperature. 
Instead, “hs intaərior motion affects the mixed layer 
temperature indirectly through the variables h, AT and We. 
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Figure 2.2 Initial Temperature Profile. 
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Pewee l veNpsra-urze-przocils is given in Figures 
202. The zeference depth, D, is 209 m and will be us=ä in 


soma of tke equations ir latsr sections. 


2. Sca Analysis and Simplification of the 


Consider now the relative magnitudes of the terms of 
equation (2.4). In the mixed lay2r, “he tenpszatu 
Change by as much as 10°C in a halt year. The time rate of 
change of mixed layer temperature would th 
of 1076 *C/sec or less. Below the nixed 4 
of change or +zmpe £ 3 ci 
be smaller. 

The first term on the right of equation (2.4) is the 
turbulent heat flux term. In the mixed layer, the only tern 
on the right side of equation (2.5) is t 


Doe heat‏ وام 
flux. Thus this term will scale with tk €‏ 


seasonal time zaze 


(I) 


of change of mixed layer temperatura, 10-6 °C/séc. The 
۳ 2 60 surface heat flux bound2ry condition, Q/ p.Cp, has 


extreme values of £4*«1073°C cm/sec. Phe value or h varies 
from around 30 m to around 100 m. Tag], tha first term en 
the right side cf equation (2.9) has a range of values 


between ‘1076 °C/sec. The second <erm on the right of equa- 
tion (2.9) has roughly the sam 


(0 


order of magnitude, but it 
can only be negative. Below the mixed layer and entrainment 
zone, chat is, below the flux region, the turbulent heat 
flux term vanishes. 

Next consider the advection “rm in equation (2.4). 
Recall that “his term is zero in the mixed layer. The 
gradient of the temperature profilə just below the mixed 
layer is initially about 1073°C/cn. During the heating 
ssason zhe temperature at the top of the thermocline will 
increase, but the vertical extent of the <hermocline will 


also increase as the mixed layer shallows. This tern should 
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thus ke expected to have the same crd=r of magnitude. The 
۶3۳۱ 1: 9۶ ٩ will be prescribed according to equation 
12:519) in section D. The maximum value at 200 n depends on 
the wave or oscillation frequancy, bu“ it will be on the 
order cf 1073 cm/sec. This ampiitads will be smaller at 
lesser depths. The result is that the advection term of 
equation (2.4) should be of the order of 10-6 *C/sec below 
the mixed layer. 

Finally, consider the diffusion term of equation 
(2.4). This term vanishes in the mixed layer. Below the 
mixed layer a reasonable scale for th2 second derivative of 
temperature is on the order of 4*C/(5000 cm) ?. typical 
range of values for the eddy diffusivity, Ss OI o 
0.5 cm?/sec. Thus, the diffusion term has a magnitude of 
about 1078 to 1077? “C/sec. The diffusion erm is much 
Smaller than the other terms of =quaticn (2.4) and will, 
from now on, be neglected. The heat equation (2.4) thus can 
ke simplified to 

dT 


— 
سے رحس —— 
ame‏ 


(2.10) 


2 
3 
1 
a 

| 


Q, 
Qu 
N 
Q. 
fu 


C. MIXED LAYER DEPTH 


In the model əf Garwood (1977), the mixed layer either 
shallicws to a retreat depth, hr, or deepens at the entrain- 
ment rate, WS. The retreat depth is determined by assuming 
steady state conditions for the turbulent kinetic energy 
budget, letting entrainment heat flux vanish and letting 
<w°<> approach zero. The resulting equation for retreat 
depth is 
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3/2 (2.11) 


dg («3/2) (T y 


(Garwcod, 1977). In this equation, m, and m, are dimension 


ma du 3 - yn {E> 
) 


0 


less model constants, « 1s the thermal expansion coefficient 
and <E> is the mixed layer mean turbulent kinetic energ 
T 


8156©, ,ينا‎ the fzicticr velocity, is the square roct 


+ O 

ct 
a 

۱-۰ 
(0 un mM 


surface wind stress divided by density. Througkou 
study, it is set equal to one centimeter per second. Sune 
E 


mixed layer mean turbulent kinetic energy Geperds on surfe 


(1) 


boundary conditions, the retreat depth is independent of the 


"O 


prescribed interior mction. 

In the case of active entrainment, the mixed layer deep- 

ening rate in the model depends on tha entrainment rate and 

Wi, the mean vertical velocity at the bottom of the mixed 

layer due to the prescribed interior motion: 
EE 7 


We - W 


= 2 
E = CPE) 


The enrrainment rate, We, derends on the entrainment zone 


(D 
۳1 


temperature jump, AT, and the 2ntrainment heat flux (see 


equation (2.8)). The entrainment heat flux is modeiled by 


1/2 ES 
mu (ww?) <E) 


(22131 
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$ 
一 | 
et 
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(Garwcod, 1977), where my is a model constant. The 
entrainment rate can then be expressed as 
AA y A DI, 
sum ( ww! E <E 
w. = pe SE (2. 14) 
xqh At 


fu 
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D com -quations (2.12) and (2.14) give a soluticr for 
dh/at if R,, <w’2> and XE» are known. 


For moderately strong wind stress, mixed layer mean 
turbulent kinetic energy, <E>, depends primarily or surface 
wind stress, which is held constant here. Ths Vertical 
component of turbulent kinetic energy, (w:2) 1-2, depends 


primarily cn surface heat flux and surface wind stress. 

When the surface heat flux becomes large enough, this term 

becomes zero, s*opping tnz entrainment process anr pernit- 

Beng shallowing. These two terms are unaffected b 
erior motion. ۱۱۳ ENN race is aii 

mean vertical motion indirectly es a result of chan 

mixed layer depth, h, and entrainment zone temperatures jump, 


AT. 


D. INTERIOR VERTICAL MOTION 


Consider for the moment a two layer, deep ocean with the 
upper layer about 200 m deep ani the lower layer about 
5000 m deep. Assuming the fluid is incompressible, the 


T 
it; 


continuity equation is difizrentiated with respect to z: 








9^ — = E = 
ما‎ = -° E EO 2 (2. 15) 
oz PETES dy dz 


where U, V and W are the msan velocity components in either 
layer. Furthermore, assume that horizontal momentum is 


vertically well mixed in the upper layer. Thus, 


۳ 
|< 


= ۱ ae: (2. 16) 
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ne 
Qu 
te 





Then W must be linear with depth in the upper layer. With a 
52055321383 appeoximaticn, Y = Oat ths surface. 

Our intersst is confined to ths upper 200 m and a linear 
vertical profile for W is a gocd approximation for long 
pericd interior motion. The functional dependence of Y on z 


is 


w(z)= - — = و‎ 2517) 


where D = 209 


m is the reference depth and Wy is the 
interior vertical motion at 200 m given by a sinusiodal 


unction of time, name 
mane tion of namel 


INE | 
W, = sin (E + ( (2. 18) 


The amplitude, W,, is defined by 








MD 15 
Wo = $ ; 2 5 23 19) 
19 lom 
For the one-year pericd interior motion, 


No = 8:10-4 cm/sec = 69 cn/day. The value of W, will be 2, 
3 and 4 times this for the one-half, one-third and one- 
quarter yeaz periods, respectively. These values of W, are 
selected so that tie displatement amplitude at 75 m is 
approximately 15 m. Ths actual values of the displacement 
amplitude depend on phase and will b3 4iscussad below. 
Ccnbinirg equations (27017) awd (25:08) yields the 


ea Kon Écr the prescribsd interior motion; 





Ww(s,t) s &- w sin( + 6). (2. 20) 
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Noites For ch= parcelle position due to ths interior 





motion can be obtained by noting that W = dz/d:, and ís 
given by 
so) ۳۱۸ ۵ ۱ iu 
z= 2, exe Le T ++) -0o:(5)) | ۰ (2220) 


For any fixed phase the maximum and minimum values of z are 


given by 
-T Vy i l 
Mar, = neo Se ( - | - Cos p) | (22727 ) 
and 1 
2 A E Nc (+ Le e) $ (2.23) 
rays DD 
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The maximum peak-to-peak displacement amplitude will cccur 
when cosg=+1 and the minimum when cos g =-1. For the 
vertical mction amplitude given by equation (2.19), a 
particle with an initial position at z=-75 m will experience 
a displacemant amplitud= (half of the peak-to-peak value) 
tetween 12.4 and 18.4 m. 

Ths model tima step is one hour. This is much smaller 
than the cne-year or even one-quarter year periods consid- 
ered. Hence, sincs Wp does not vary significantly over a 
cre hour time step, “he integral of (2.17) can be approxi- 


mated by 





at) (2525) 


where At=t,-t, is the time increment (one hour) and VW», is 


9003532526 4 from equatien (2.18). In the model, vertical 
advection is handled by compressing or expanding the 
Mert acal grid according to equation (2.25). No significant 


difference in tha cesuits was found between model runs using 
equation (2.25) and equation (2.21). 

The vertical increment is initially one meter, but this 
varies as aresult of the compression or expansion of the 
Meeeical grid. This variation in the vertical increment is 
never more “han about 25%. A consequence of this scheme is 
that the lowest grid poin- is not always a- a depth of 199 m 


Ben the initial prefile, Figure 2.2. Care must be taken, 
therefore, when performing vertical integration from the 
surface to the reference depth, D = 200 m. The integrals 
that will te calculated in cha nex z section invclve “he 
temperature differance, T-T, - Probiems with varying the 
grid size can easily be avoi idea by setting che reference 
temperature, T,, qual to Tp, the initial temperature at 
200 m. This way, T-T, will always be zero at the base of 
the vertical grid. The grid is initially made deep enough 


so that the seasonal thermocline never reaches the bottom of 
the grid; hence, T, always equals its initial value, 6.08°C. 
Subtracting the reference temperature, T,, from tempera- 
ture, T, in the integrals in tha naxt section amounts to 
Subtracting the terms: 
0 


EN T D (2. 26) 
~D 


and 
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= = T و‎ (2.27) 





] + ان ج‎ 2 
D) 
=") As 
Since T, always remains constant, Sth expressions, 2,26) 
end (2.27), will always remain constant. Sera coa T. 


from T will amount to subtracting constants from the 
integrals ir the next secticn. 


E. POTENTIAL ENERGY AND HEAT 


The potential energy per unit area for a column of water 


from the surface to depth, D, is given by 





PE £ | (2.28) 
-). gpa dez. 


(LI E Co 


ia 


Employing equation (2.2) with T, as the reference 


temperature, 2llows equation (2.28) to be written as 


١ TAE 2 0 | | TOUR 
Et Po AC ec ( PE i) 2 AE + Constant. ۰ 








e 


=) 


Thé constant in equation (2.29) is 








| 4p, 2 da EA De: (2. 30) 


ES Po X. (X 


The integral in equation (2.29) will be used for the compar- 
[on Om the potential energy with and without the interior 


motion. It wili be represented by the variable P with units 


ore me wher 


P= -Í ۱ A " (2:31) 
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A RSA: content vilt ba indicated by the variable H 


with units °C m whers 
o. (2232) 
En a+. 


The analysis of the physical effects of the interior 
motion can te facilitated by developing the expressions for 
the time rate of change of H anid of P. ۳۰ 


differentiating equation (2.32) with respect to time gives 


9d 3 | (=. =" DE) 
DEC E! (3 >) de 


Equation (2.33) can be written as 


a (2.34) 


OH ۳ o 


N 
Due 


y 
c t - D 


Substituting equation (2.10) for th= time rate of change of 


temperature into equation (2.34) yields 


H i y as M 9 
له‎ A fea em 
yc En = 25 c 





Q_ 


The first term on the right, when integrated, is just the 
het surface heat flux. Spe cleucang equa=ion (2.17) ¡nto 
the second term on the right of (2.35) and integrating by 
parts gives 
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(2.535) 




















m سس‎ y = Ww 7G 
E Wide -D[-vt > | Tar) : 
-D 97 ! -ò 
Using equation (2.26) this can be réwritten as 
f E J TE (23 37) 
- Ww — da = MD = 一 个 Lom .; 
ھت ول‎ Bra = en 
-D 5 
۲ ۱ 0:5 (2.35) can new be written as 
EM —G - We y (2. 38) 
E Po Ey D 
Me” tommulation for the time rate of change of the 
potential energy is Similar and is given by 
a حم‎ C y e 
el = - { c> i 2 le: (2.39) 
JE ct 
-D 
5۱95 ۰2۱15: 280 equation (2.10) into (2.39) yields 
s «2 ` c \= 
SER + 29 (wr) eds + | 2 W EE. A? (2.40) 
J c) à E 
M -D -d ١ 


in equazion (2.40) can be integrated by 


115 first term 
parts: 
۰0۵ si a م‎ 
|: SG) ds (EF) Å ET) aa. (2.41) 
له‎ 
e م‎ 


D 
of squation (2.41) vanishes 


W'T' is zero peicw the 


term or the right 


= 0 at z= =D. 


The forst 
because w'T' 


Since 
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entrainment zone, equation (2.41) reduces to 


EAD. 1 m (2002) 

| E le — | (wT) dz. 

et - 
Rc 


This integral may ba broken up over the mixed layer and the 


entrainment zone, prcducing 


C = 


P mh... c (2408) 
۱ 2 MTT) د‎ wi Tl A® و‎ 
Eo: ۱ 
“hi 
wnsre <w‘T’> is the vertically averaged mixed layer mean of 


the turbulent heat flux. Since w‘'T’ is linear with respect 
to z ir the mixed layer, this mixed layer mean can be 
written 





-aT st + AT We). (2.44) 


Consider the seccnd term on tae right side of equation 
(2240) . Substituting equation (2.17) for W and integrating 
by parts yields 





using equation (2.27) this can be rewritten as 


> < چاه له | 


o) 
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O 
ل‎ Dm 
۱ >| Ge) ee SAA (2. 46) 
> 2 D 


On combining 2quations (2.40), (2.43), (2.44) and (2.46), 


the time rate of change of the potsnzial energy variable now 
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becomes 


r | 
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the turbulent EX 


entrainment zone, as illustrated in Figure 2. 1: 
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Because is monotonic across 


ch 
| ST NIT Pee A Wwe 


-h-ó 

This indicates char the (2 
Should genscally bpe smaller in magnitude than the first 
the right 


entrainment 


integral tarm in equation 


en side, unless of course, tne surface 


heat fluxes nearly cancel. 


Actual values of the terms of C103 (2.47) 


Calculated for one-year motion ac 


different phases and ths the right of 


equation was smaller than 
magnitude. Furthermore, 
generally an order of magni-ude 


other two. Equation (2.47) will th 
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III. EXPERIMENTS AND RESULTS 


A. SUMMARY OF PERTINENT EQUATIONS 

















The =quations for time rate of change of mixed layer 
depth, mixed layer (and sea surfacs) temperature, heat and 
potential energy are 

sh GE (3.1) 
C , ‘ ۰ 
GL G AY AT We (3.2) 
«we تست .( جر‎ —————— ? 
E Po Cp h h 
v H G LV | 
mo = - 2H و‎ (3.3) 
EE Po cp D . > 
and 
۱ e | 5 د‎ (3.4) 
ELO i , 
Et 2 PCp 2 D 
The retreat deptk is given by 
l 2 . — mea 
h, = Ma y + m, SE) (3 .5) 
AS. ss. ۰ 
a o "n 
= p 


where (W'T'), has been replaced by equation (2.6). 


The entrainment rate is given by 
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The prescribed net surface heat flux is given by 





G e n EN \ (Sr 79 
Po Cp Pop | year 


and is plotted as a function of time in Figu 


Q/ (> 2 ESE Geese) 
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Figure 3.1 Net Surface Heat Flux. 


The prescribed vertical motion at the reference depth, 


REC (3.8) 
W > = W. Sin [5D + +) i 
where W, is defined Ly equation (2.19). By Substituting 
z=-h into equation (2.17), it can be seen that W, varies 


Gamectiy as h and W,. 
Some additional nctation will be useful for the discus- 


sion of the results. To observe tha effect of the interior 
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۶۶ ۰18۰ the values of he Tus. -H aad > are compared to those 
enouz interior mction. conteos ase plotted ct “ha 
difference of mixed layer depta with and without the 
interior motion, likewise for mixed layer temperature, heat 
and potential energy. Therefore, the following new 


Variakles are defined: 


A |, —- h | (3.9) 
(3. 10) 
O 


AP EDO E (2. 12) 


B. RESULTS WITH NO INTERIOR MOTION 


The mixed layer depth has been plotted in Figure A.4. 
Note that the vertical axis scale is in 
The year was started at the beginning of the heating season 
and the mixed layer retreats for the f 
Then the mixed layer deepens slowly for the 
and mor2 rapidly in the last quarter. In the last few days 
the mixed layer shallows tc the <=0,1 year,...,n-years 


cyclical steady state value of about 75 m. 
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The mixed layer temperature, T is plotted in Fig 


a ur 
ASD + As expected, the mixed layer temperature inc 
during the heating seascn and decrsases during the cocling 
season, for the most part. However, ncts that th 

mixed layer temperature occurs around day 170 and T,, starts 
to decrease ten days befere ths ens cf the heating season 


9 


This is because the entrainment hea 


ct 


flux term has begun to 
dominate equatio v2). Some entrainment had already 
Started fifty days earlier. By day 170, towards the end of 
the heating season, AT has reached its largest value, 
NES C. Cn day 170 the cooling of the mixed layer due to 
entrainment has reached 0. 021 °C per day while che surface 
heating has dropped to 0.030°C per day. 

The heat content, plotted in Figure A.6, depends only on 
surface heat flux for *he casa of no interior motion (see 
equation (3.3)). Any heat lost from che mixed layer dus to 
entrainment will be gained by the region below. Since hea+ 
is calculated from the surface to 200 m, entrainment will 
have no effect on its value. 

Potential energy per unit area is plotted in Figure A.7. 
This, “oo, generally increases with heating and decreases 
sach cooling. However, note that the peak value occurs at 
about day 305, 25 days into the cooling season. This is 
because potential energy is also increased by mixing as 
indicated by equation (3.4). The entrainment term dominates 
equation (3.4) from day 170 until day 305, wnen the net 
surface keat flux is near zero. A chose imepection. of 
Figure A.7 reveals that the potencial energy at the end of 
the year is greater than the initial value. Spec: fical lys 
the initial value of E is 19608*C m?; che value at the end 
of the first model year is 20204°C m? and at the end of che 
second year is 20579°C m2. This increase in the potential 
energy is due to the adjustment? cf the temperature profile 


as the ocean approaches cyclical steady stats. 
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Figure 4.1 is a plot cf tempsrature versus depth and 


tine. The top cf tte graph where the temperature contcurs 
are vertical is che mixed layer. The mixed layer structure 


is quite evider- after the first quarter year due to the 


0 
strong temperature gradient at the base of the mixed layer. 
Note that the temperature profile at any depth below the 
mixed layer is left unchanged until mixing reaches that 


depth. 


C. ONE-YEAR PERIOD INTERIOR MOTION 


The zero contours of ha prescribed interior motion have 
keen drawn in Figure A.8 with upwelling and downwelling 
indicated. The dotted contours indicate when the maqnitude 
of the interior moticn drops to half of the maximum value. 

During “ha first quarter year the mixed laver is skal- 
lowing. Since the retreat depth depends primarily or 
surface Ecundary conditions, mixed layer depth will be unaf- 
fected by interior motion at this time. This is borne out 
eye ome contours of Ah (defined in equation (3.9)), Figure 
MC. The first quarter year shows no difference in mixed 
layer depth. With little cr no effect on mixed I 

will 
difference in the value of the first term 


layer depth 
during “he first quarter year, there bp liete "cro 
of equation (3.2) 
WEENM end without the intericr motion. The  ssccni term of 
equation (3.2) will te zero at this time Since entrainment 
Stops during shallowing. Hence, mixed layer temperature 
will alse remain unaffected by the interior motion during 
the first quarter year. This can be seen in Figure A. 12. 
Heat and potential energy, on the other band, are influenced 
by the advection of temperatura in the <hermocline as indi- 
cated by she presence of the last term in equation (3.3) and 
equation (3.4). Thus, he intemior Motion will influence 
hear corSenfeard potential energy even during the first 


quarter year (Figures 4.16 and A.13). 
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1. Effect on Mixed Layer Depzh 
evar. On SCS ..3) indicates that downwsllina will 
enhance deepening, and upwelling will reduce deepening. T 


Figure A.9, i< can be seen that dsspening still ces occur 
after the first quarter year as iz did withcut any interior 
motion. However, the rate of deepening depends on phase. 

At the 90° phase, there is upwelling during the 
middle half of the mcdel y2ar (ses Figure A.2) reducing «he 
rate of Geepening. At day 274 downwslling begins and the 
mixed layer bagins tc despen more rapidly. The very rapid 
deepening in the last quarter of tas year for This chase 
results Írcm both antrairment and downwelling working 
together. Calculated vaiues of the terms of equation (3.1) 
indicate that entrainment is the dominant term most of the 
time during deepening. In this example zhe entrainment rate 
has been enhanced by conditions sst-up by the interior 
Moticr. Equation (3.6) indicates that entrainment rate, We, 


is inversely proportionai to mixed layer depth, h and 


= 


entrainment zone temperature jump, AT. Compare Fiqure A.2 
COM QuUre امش‎ (Ste case with no interior motior) and A. 3 
(the 270° hass case). It can be séan that at the start of 


the las: quarter of the year the 90° phase case has a much 
ur 


(D 
ct 
3 
[n 
3 
(1 
SF 
iD 
0 


lower mixed layer remperat zher two cases (the 
1 


reason for this will be discussed in the next section). 


This lower mixed layer temperatur2 means that the entrain- 


ment zone temperature jump wili be smailer, allowing for a 
larger entrainment rate. At the start of the last quarter 


of the year, the differences in h are not as significant as 
those in AT. 


For example, the values of AT, h, We and Y, have 


teen calculated for these three cases for dav 280. The 
mixed layer depth, he is about 50 m for the case without 
ineerio- motion, 65 m for the 90° phase and 68 m for the 
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270? phase. These valuss de not greatly differ. However, 
the temperature jump for the 90° phass case is 1.3°C which 
35 bess than half the 2.9°C value for the case with no 

than the 4.7°C value 


for the 270° case. The entrainment rate for the 90° phase 


h 
9 
intericr motion which in turn is less 
2 
to 43 cn/3 


case is 90 cw/day compar 2d a 
motion case and 22 cm/day for tha 27 


y for -he 20 inserior 
0° case. For "Boch 
phases considered, the interior motion is at 1 
of magnituds smaller than tn 
Neverthslsss, at day 280, the d=2psning in ti 0 

case is helped a little by downweliing and in the 270° case, 
hindered a little by upweiling. 

In the last few days of the year, “he increase in 
mixed layer depth begins co take over, 
90? cass. This inhibits further antca 
layer then rapidly shallcws to the cyc 
value of 75 m (Figure A.9). 


In the 270° vhass cass (Figure A.3) there is down- 
weliing during the middle half of <he year, which increases 
the deepening of the mixed layer. By the start of the last 


quarter of the year, “he mixed layer is warmer than in the 
case with no interior motion or in ths 90° phase case. This 
leads to a larger value of AT and a reduced entrainn 


rats. 


Hb 


265-3159 Figures A.8 and A.9 it can be seen that 
mid-year upwelling is associated with 2 20 

deepening and mid-year downwaliing is 
increased deepening rate. Looking at Figure A.10, it can be 
seen that the actual decreased or increased valu 

layer depth lag the start of upwelling or downwelling 
(respectively) by 90°. The time rate of change of mixed 
layer depth tracks alon wii ERS interior motion in 
accordance with equation (3.1). However, actual mixed layer 


0۳5 Ch the integral of this and integrals of periodic 
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functions lag them by 90°. Also note that ths strong 
enhancement of and strong hindrance to deepening at “ha end 
of the year (Figure A.10) lags the strong decrease and 
strong increase in mixed layer temperature (Figure A.12) by 
about 90°. 


te pefe On Gaxsd Layer Temperatura 


Consider now the mixed layer heat equati 


0 
The net surface heat flux is prescrined and unaffected by 


as interior notion. The entrainmsnr. neat flux, giver by 
equation (3.6), is inversely proportional <o mixed layer 


depth. The other terms in equation (2.6 


) 
sem So, che fir 


on surface boundary condition - 2 120 2 sequa- 
tion (3.2) is inversely proportionai to mixed layer depth, 
h, and the second, inversely proportional to the square of 
ha This implies that ar increase in h will reduce “hs 
Magnitude of mixed layer temperature change; a décrease in h 
will increase the change in Ty. 1۳ وک‎ became, wishes 
deeper mixed layer, a given amount of heat added + QE 
d 


removed from tha mixed layer must ba sprea 
mass cf water. 

ct is illustrated in “ns examples shown in 
Figures A.2 and À.3. 1 0 ۲ ع ح ناى‎ A.2, the mex 


ed 

shallower than in Figure A. 3 from about day 190 to abou- day 
E 
th 


300. This resulted from the mid-y22r upwelling for the 90° 
phase case and the aid-year downwelling for the 270° casa. 
The shallower mixed layer in Figure A.2 allows for a much 
more rarid temperature drop than occurs with tne deeper 
mixed layer of Figure 4.3. This can bs seen in the enhanced 


h = 
gradients in Figure 4.2. After day 300, che rapid deepening 
of the mixed layer for the 90° phase case andi slow deep 
for the 270° case reverse the situation. 

To see the effect of the interior motion 


f 
the phases, the Ta contours should no: be compared dire 
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ET OG 5502 eantcurs but rather to the contours 
of the mixed layer depth difference, Ah (defined by equa- 
tion (3.9)). Between the 0° and 189° phases, after day 180 
these is a region in Figure 34.10 where the mixed layer is 
shallower with the interior motion than without. mhis 
causes the mixed layer temperature decrease (the second half 
of the year is the ccoling season) to be greater, resulting 
in the negative values of AT,, in the last quarter year at 
these phases (Figure A.12). Between che 180° and 360° 
phases the opposite occurs. 

At the end of the year the mixed layer depth returns 
ED 75 M. Mixed layer temperatura does not return to its 
initial value. The second year T,, data is plotted in 
Figures A.13 and 8. The reduced temperatures for phases 
0° +o 180° anã increased temperatures for phases 180° to 
360° are continued nto the second year The mixed layer 


depth follows the same pattern in the second year as i+ did 


En the first. This causes the T,, pattern begun in the 
first year to be enhanced in the sacond year. Continued 


pumping of “ha ocean by this period interior motion enhances 
this structures. If permitted, i+ would eventually cause the 
upper 200 m to become isothermal (in the region of the 90° 
phase ar the end of the yaar), followed by deepening beyond 
the limits cf what the model can handle. It would not be 
realistic, however, to have this single pario interior 
motion continue for the several y2ars necessary for this +o 
happen. ee <The interior notior would probably be 


damped by interactions witn “he boundary layer mixing 


processes. There should 4 be more than one period 
present. Furthermore, surface heat flux should be infiu- 


enced Ly sea surface temperature. These three possibilities 
se 


open up avenues for further 
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3. Eiíiscz on Heat Content 
Equation (3.3) indicates that change in the heat 


content depends on net surface heat flux and on pimping cf 
the ocean by tae interior motion. Ths net surface heat flux 
is prescribed and is independenz of the inter: 





9 
Wib-racting equation (3.3) TOF ER: casa Of RO interior 
motion from the same equation with interior motion qives 
ETT, (3. 13) 
ci © D 


where AH is defined by equation (3.11). Egyar.on ( 
indicates that upwelling will cause AH to decrease and 
downwelling will cause AB “o increase. Those 
cbserved by comparing Figures A.8 and A.16. The input to 
heat content due to net surface haat flux alone causes H to 
increase during the first half of the year and dec 
during the second half as shown in Figure A.6. The com 
tion of surface heating with tha affect of the int 
motion produces the saddle-shaped curve shcwn in Figure 
AND, 

The effect of different phases of the i 
motion can ke demonstrated by the following two simple ex 
Més (Figures 3.2 and 3.3). No mixing is permitted an 
two-layer ocean is assumed for simpilci <y. 

In “he first example shown in Figure 3.2 thera is 


initial downwelling followed by heating, cher upweliing 
E Sa 


Co eme Ey cooling and final downwelling =5 tha initíal 
mixed layer davth. when do — os precedes p the 
heat is distributed cver a thicker mixed layer. A fixed 


amount of heat wili produce a smaller increase in T than 
wcuid have resulted without the downw2lling. DL ^ tho ds 


followed by upweiling before cooling, chè thinner mixed 
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Figure 3.2 Heat Content Example 1. 
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Figure 3.3 Heat Content Example 2. 


43 





layer will reach a lower temperature for a given amount ct 
heat removed. At the end of the complete cycle the mixed 
layer temperature will be lower. Recall that the hear 


content is defined by 


hē lower mixed Layer temperature produces a lover 


Upwelling preceding the heating season distributes 
C 


the heat cver a “hinner mixed layer. This allows a greater 
increase inT than would otherwise occur. Tr this is 
followed by downwelling preceding the cooling season, a 
fixed amount of heat removed wil produce a smaller drop in 


T. The net effect will be a higher temperature in the mixed 
layer and a higher heat content. 

Because of the asymmetry, as illustrated in these 
examples, heat is nct conserved for = particular phase of 
۳ هد‎ (02 0 This can be seen in Figures A.15 and 
2.16 . at some phases the year-eni heat content is higher 
than it was at the beginning of th2 y2ar, at other phases 
lower. However, if the heat content is averaged over all of 
the phases, the year-end vaiue equals the initial value. 
The average heat content for sixteen equally spaced phases 
from 0° to 360° was tabulated. The phase-averaged value of 
H was 0.07°C n less at the end of the year than it was 
initially. This small error can bs aviributeã +o the coarse 

1 


resoluticn of averaging only sixteen phases over the entire 


360°. 
4. Effect on Potential Eneray 
Time rate of change of potantial energy is given by 
equation (3.4). From this equation it can be seen that 


Surface heating *ends to increase — پوو‎ energy, surface 


44 





> Au 








TA ۰ 
۱ Å iT tow 


cooling tends =o decrease it. The net surface heat flux is 
independent of the interior motion so the effect on ¡AP due 
to this *erm depends on dah. If there is an increase in 
mixed layer depth during heating, the buovancy added wili be 
mixed deeper, increasing potential energy. If there is an 
increase in mixed layer depth during cooling the effect will 
be opposite. 


The entrainment heat flux term has a positive sign 


in equation (3.4). Increased mixing moves buoyant water 
down and less buoyant water upward, increasing potential 
energy. 


Ccmparison of values of the rerms of equation (3.4) 
indicates hat the third term on the right is generally 
dominant. The exception is when the interior moticn goes to 
zero, twice a year. The first term of equation (3.4) will 
not differ greatiy with and wichcut interior motion unless 
there is a larga difference in mixed layer depth. Recall 
that entrainment heat flux is nearly inversely proporticnal 
to mixed layer depth (equation (3.6)). Thus, the second 
SEM will nct vary significantly. SO, the potential energy 
closely follows the interior motion. This can be seer by. 
comparing Figures A.8 and Aa 18. Durirg dowmmelling, 
potential energy increases, during upwelling i= decreases 


The physical reasoning behind the effect cf “he 


third term of equation (3.4), the pumping tern, is that 
downwelling pushes buoyant surface water downward, whereas, 
upwelling is in the same direction as the buoyant forces. 


2. Semma-y otf esponse to One-Year Period Interior 


During the first quarter sar there is shallowing 


Y 
ad the 1nterior:fotion hás little effect 
å 


on either mixed 
layer depth or temperature. Mid-year upwelling reduces 
deepening and mid-year downwelling enhances deepening. The 
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shallow=r mixed layer MES cuina from upwelling wil eco! 
faster; he desper mixed layer resulting from downwelling 
will cool slower. A cooler mixed layer has & reduced AT 
allcwing for a larger entrainment rate and more rapid 
deepening. A warmer mixed layer despens more slowly. 

Potential enezay and heat are both directly affected 
by the interior moticn. Beth are increased by downwelling 
and decreased by upwelling. 

Before examining the mixing response tc higher 
frequencies, an observation can be nade regardira hcrizental 
eezucture produced by the interior motion. As noted in 
chapter one, phase may be related to horizontal displacement 


by substituting 


)14 .3( ._ ر 


into equation (2.20) to obtain 








we: +: W. sen (ATL + 1۲ ۰ Gane 


The wavslenath, L, might be on the order of 300 km, large 
enough for the one-dimensional assumption to be valid. 
Figures A.12 and A. 14 show ths development of hori- 
zontal bands of positive and negaziv= mixed layer tempera- 
ture anomalies. which would result from a one-year period 
interior wave. Because temperatura is assumed to be verti- 
cally homogeneous in the mixed layer, T,, also represents 
sea surface temperature. Figures A. 11 and A.13 show a time 


sequence of the horizontal sea surface temperature field 


nscuTting from" the interior motion, the heat cycie and 
mixing processes. The entire surface warms during the 
heating season and cools during the cooling season. Dés ing 
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MESES tt RNA 62 the First year, there is no significan: 
difference in 


» 

ea surface temperature along the horizontal 
year, the region where there is 
È 


5 
axis. Aces ER2 first half 
upwelling at the start of the cooling season becomes and 
remains cooler than the region where there is downwelling at 


the start cf the cooling season. 


D. EFFECT OF SHORTER PERIOD INTERIOR MOTION 


Ccmparison of Figures A.19 and A.21 shows that as before 
downwelling increases 62 202 23 5 and upwelling reduces deep- 
ening. This can also be observed by comparing Figures A.30 
E32 ana by comparing Figures ۸۰ 39 and A.41. The more 
rapid variation in mixed layer depth does not yield the 
strong effect on temperature found for the one-year period 
۳۱ 3 ۵۲ motion. This can be seen in the weaker gradients 
of AT, cf Figures A.23, A.34 and A.43 compared to that of 
nM. Hence, the higher frequency interior moticn does not 
hav= as strong an affect on the mixed layer as the lower 
frequency mcticn does. 

The second year temperatuzs data was plotted for the 
half-year period interior motion in Figures A.24 and A.25. 
y 


å 


ct 


Notice chat the same structure sezn in the fi 


1 


c 
ta 


(0 


ar reav- 


(D 
E 


pears in the second, but the cold spots hav ot grcwn 
2 


significantly colder nor the warm spots warmer a in the 
case oft tke one-yeaz pericd motior. 

As before, the heat and potential energy differences, 
AH and AP, track along with the in-erior motion as seen in 
۳۶ ۷ 2 65 A.27, 25292 “Beso, BE, M5 and ww q7. The 
resulting contours of H and P apvear more complex. This is 
because *he values cf H and P without the interior motion 
follow a one year cycle (Figures A.5 and A.7) and the higher 
frequency cycles due to the pumping by “the intericr metion 
are superimpcsed on the one year cycle. 
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Fcr the cne-year period interior motion, upwelling lasts 
long enough to enhance cooling significantly (at phases 
around 90°) and downwelling lasts long enough to inhibit 
cooling significantly (at phases around 270°). Themes ne. ter 
periods de not maintain upwelling and downwelling long 


enough o produces an effect this strong. 


E. ASYMMETRY IN ENHANCED AND REDUCED MIXING 


In Figures A.48 and A.49 the values of peak enhanced 
cooling and peak reduced mixed layer cooling are plotted. 


Both first and second year values ars considered (in Figures 
A.48 and A.49, respectively). Negative values of AT. 
indicate enhanced cocling due to “he interior motion. The 
maximum abeclute value of all of tht negative Afu, data 
points was determined to indicate the peak enhanced cooling. 
Positive values of AT. indicata reduced cooling. The 
Maximum value of all of the positive AT,, data points indi- 
cates the peak reduced cooling. The positive valu2s do not 
generally indicate enhanced heating, because the wav2 has 
very little effect on the mixed layer temperature during the 
heating season. 

The first obvious observation that can be made from 
Figures A.48 anã A.49 is that the largest peak anomalies 
occur with the an=s-ysar period interior motion. For the 
cne-year period interior motion the paak enhanced and peak 
reduced mixed layer cooling is 2.04 and 2.16%, -=spec- 
tively. For the one-half year period motion *hese values 
are 1.35 and 1.15, respectively. Meco” “for each of the 
periods studied, the difference between the peak enhanced 
ccoling andthe peak reduced ccoling is greater in che 
second year than in the first. All of the periods studied 
except tte one-year period hav2 greater peak enhanced 
cooling than peak reduced cooling. 
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The phass-avsraged values of the variables ah, AT, 


= 


AH, and AP were computed for sixtesn equally spaced phases 


EM 0° to 360°. The phass-averagsd values at the end of 
the first and the end of the second model years are listed 
Me. These values indicate the nex effect of the 
interior motion after e2ach year of interacting with the 
maxing processes. Consistent with the results shown in 


Figures A.48 and A.49, the results in Table I indicate that 


bi 


the one-year period interior motion has produced a net 
ređuctior of mixed layer cooling. The higher frequency 
interior motion has produced a net enhancement of mixed 
layer cooling. For all of the pariods studied there has 
been a net reduction of the mixed layer depth. Heat is 
conserved, hence the phase-everagsd value of AH is zero at 
the end of each year. The effect of the one-year period 
interior motion was to decrease the potential energy. The 
higher frequency interior motion produced a net increase in 
potential energy. 

A net increase (decrease) in potential 2nergy indicates 
that the interior motion has enhanced (reduced) the overall 
me xing. The amount of deepening required to produce an 
equivalent increase in potential energy is instructive. 
This can easily be calculated for the simplified example of 
a two-layer ocean, illustrated in Figure 3.4 with initial 
mixed-layer temperature, T, and depth, h,. This mixed layer 
is deepened to h, with new mixed layer temperature, T.. 
Heat is censerved, so, using the definition for tha heat 


variakle, H, given in equation (2.32): 
or 3 Jh, = ne 1 Ie ) hs سد‎ 


1 


Using equation (2.31), ¿ns change in the potential 


m 
"7 
(D 
ti 

2 

ng 
pa 
n 


given by 
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Substituting 2quation (3.16) 


Dal 
Za 


CS 5223 T,-T, is just the 
The 


heat 


ture jung. increa se 


ening with conserved 


increase in 






initial profils 


Figure 3.4 


Consider 2 
a+ the 


the 


lies of the 


Suppose 


intial re jump 


potential energy for 


would then ke 


mixed layer depth 


the phase-averaged 
fis 
initial gi xed 
u ism °C, 
the one-year 


equivalent tea 


(90%) 


into equation (3.17) yields 


AQ ( (3. 18) 


initial entrainment zone tenpeza- 
in potential energy due to deep- 
to 


2xample. 


i. di =cs ly proportional 


the 


for the two-layer 


= = د = ےد 
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after mixing 


Two-Layer Ocean Example. 


potential energy anoma- 
in Table TI. 
is 756m the 


Gift 


year, shown 


layer depth and 


The net reduction 
period interior motion 


reduced mixed layer depth of 
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0.73 m. The net increase of votential energy for the on2- 


As 
Í 3 
3 
3 
(D 
tt 


half, one-third, and one-quarter yaar perio a 
is equivalent to an increase of mixed layer d=pth 
1.8 and 1.8 m, respectively. 

The amount of dcwnweiling and upwelling are the sane 
over the whole year and for ali phases. Yet, the net reduc- 
tion of mixing and mixed layer cooling does not equal +he 
net enhancemen- cf the same. The net effect of the one-year 
rodd Interior motion is to hinder mixing. The net effect 


of the shcrter period interior motion is to enhance nmixirg. 


TABLE I 


Phase-Averaged Mean Anomalies 


Phase-averaged mean of: 


Day Qh Ade AH AP Period 
(m) (°C) (^C m) (^C m) (years) 





0. 0. 
365 
739 





365 
730 


122 


365 
730 


1/3 


365 
730 


1/4 
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IV. CONCLUSTONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


Seasonai-scale interior motion has a significant effect 
on boundary layer mixing processes. Th 6 

pronounced when the period of the interior motion is cne 
year. For this pericd, upwelling between the heating and 
cooling seasons can strongly enharce mixed layer cooling and 


subsequent deepening. Downwelling can strongly reduce the 


same. This effect may produce a significant observable 
horizontal ra-tszn in sea surface temperature Although 
there is negligible effict on sea surface temperature and 


mixed layer depth during the beginning of che heating 
season, the interior motion does affect the potential energy 
and heat content at this time. Threughout the year the 


potential energy and heat are directly related to the 


pumping by “ha interior motion. Upweliing decreases and 
downwellin increases both potential energy and heat 
content. 


Higher frequency intericr motion does not have as strong 
an effect. Ths cycling betwasn upwelling and downwellina 
during the cooling season prevents the strongly enhanced 
mixed layer cooling from occurring. 

Tne effect of the interior motion is asymmetric. 


Generally, the enhancement of mixing and mixed layer ccoling 


exceeds the reduction f the same. Piles) @XGeption to this 
۳ ۱5 Cor the ona-y=ar pariod intarior motion. 
B. RECOMMENDATIONS 

One aspect of this research which departs from reality 


ES 
O 
ct 
{ dle 
O 
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ncn Seeeseribing a singl wave cor the interior 
This could be changed quit? easily by prescribing 


D 








| Ait (4.1 
Wy = بو‎ ..« ) E ro), ) 


fhe amplitudes, W,, and periods, ,يلا‎ could be based on the 
first few dcminant frequency components observed in actual 
wave spectra. 

Another way that reality Cola De more closely 
approached is by allowing tha net surface heat flux to vary 
depending cnthe sea surface temperature. This should 
dampen the strongly enhanced or reduced cooling of the mixed 
layer temperature seen with the one-ysar period motion. 

An area open to further research is to use a finite 
depth ocean and an appropriate vertical wave structure. An 
appropriate wave structure for planstary waves is given by 
Willmott and Mysak (1980). In their work the vertical 


motion is linear in the mixed layer: 


2 
LI و‎ ae 
RE 


but, below the mixed layer the vertical motion is given by 





ER y + 
© Ww 5 - 3 Wi - 0 (4.3) 
ae cos 


where N2 is the Brunt-Váisalá frequency and Cc, represents 
the speed of propagating long waves in che non-rotating cass 
for mede number n. The solutions oí equations (4.2) and 
(4.3) must be matchad at the base of *he mixed layer,  z--h, 
to ensure continuity cf the vertical velecity. The boundary 
conditions are rigid lid and no normal flow through the 
ocean bottom, so, w=0 at z=0 and at the ocean bottom. At 
each time step the mixed layer depth is recalculated and so 
the vertical wave structure problem would also have to be 


solved for the new w(Z). Tha combination of this sort of 





vertical wave structure with the Garwood model for boundary 
layer mixing processes could provids an interesting area for 
future study. 

In “his study ute interior motion was prescribed and 
acted on the mixed layer. The interior motion was not 
allowed to te changed by the mixed layer processes. if, for 
example, zune wave motion enhances mixing, the increase in 
potential energy should be balanced by a decrease in wave 
energy. Although the frequencies considered here are much 
lower than those in the work by Bell (1978), the processes 


involved could bea sirilar. 
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Figure A. 1 Temperature, No Interior Motion. 
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Figure A.2 Teaperature, One-Year Period, 90° Phase. 
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Figure A.3 Temperature, One-Yaar Period, 270° Phase. 
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Figure A.U Mixed Layer Depth, No Interior Motion. 
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Figure A.5 Mixed Layer Temperature, No Interior Moticn. 
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Figure A.6 Heat, No Interior Motion. 
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Figure A.8 One-Year Period Interior Motion. 


59 





PHASE DEG 
90 135 180 225 270 315 360 


45 





Q0 30 60 90 120,150 180 210 240 270 300 330 360 
TINE (ATS 


Figure A.9 h (meters), One-Year Period Interior Motion. 
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Figure A.10 Oh (meters), One-Year Period Interior Motion. 
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Figure ۸۰ 1 Ta (C), One-Year Period Interior Motion. 
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Figure ۸۰13 T,, (C), One-Year Period, Second Year. 
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Figure A. 4 ATu, (C), One-Year Period, Second Year. 
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Figure A.15 H (100°C m), One-Year Period Interior Motion. 
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Figure A. 16 AH (100*C m), One-Year Period Interior Motion. 
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P (10**C m2), One-Year Period Interior Motion. 
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AP (10*°C m2), One-Year Period Interior Motion. 
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Figure A. 19 1/2-Year Period Interior Motion. 
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Figure A. 22 Ta. (GR 1/2-Year Period Interior Motion. 
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(C) , 1/2-Year Period Interior Motion. 
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Figure ۸۰24 Th, (C), 1/2-Year Period, Second Year. 
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Figure A.25 AT. (©). 1/2-Year Period, Second Year. 
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Figure A.26 H )100 m), 1/2-Year Period Interior Motion. 
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Figure A.27 AH (100°C a), 1/2-Ysar Period Interior Motion. 
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A.29 AP (10% m2), 1/2-Year Period Interior Motion. 


70 





Pt ty 


۳2۱۵ 


Figure A. 30 


. T N 
È T 












ES 
a 
۵ 


> 


\ 
lr m D 
es 
TIME DAYS 


1/3-Year Period Interior Motion. 
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Figure A.31 h (meters), 1/3-Year Period Interior Motion. 
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A.32 Ah (meters), 1/3-Year Pariod Interior Motion. 
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Figure A.33 Ta, (C), 1/3-Year Period Interior Mction. 
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Figure A.34 OT. (C), 1/3-Year Period Interior Motion. 
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Figure A.35 H (100°C a), 1/3-Year Period Interior Motion. 
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Figure A.37 م‎ (104°C m2), 1/3-Year Period Interior Motion. 
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Figure A.38 AP (10%*C m2), 1/3-Yaar Period Interior Motion. 
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Figure A.40 h (meters), 1/4-Year Period Interior Motion. 
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Figure A.41 Ah (meters), 1/4-Year Period Interior Motion. 
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Figure A. 42 
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Figure A.44 H (100°C m), 1/4-Year Period Interior Motion. 
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Figure A.45 AH (100°C m), 1/4-Yaar Period Interior Motion. 
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Figure A.46 P (104°C m2), 1/4-Year Period Interior Motion. 
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Figure A.47 AP (10%° m2), 1/4-Ysar Period Interior Motion. 
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Figure A.48 Peak Enhanced and Reduced Cooling Values, Year 1. 
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Figure A.49 Peak Enhanced and Reduced Cooling Values, Year 2. 
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